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ABSTRACT 


Define  F  <  G  (F  <  G)  if  F  and  G  have  the  same 
median,  say  the  origin  and  _G  *F(x)  is  concave- 
convex  about  the  origin  (G  *F(x)/x  is  increasing 
(decreasing)  in  x  positive  (negative)).  Conservative 
tolerance  limits  are  derived  for  distributions  which 
are  s-ordered  with  respect  to  the  Laplace  distribution. 
These  are  especially  reasonable  for  mensuration  data. 

In  addition,  many  inequalities  concerning  combinations 
of  order  statistics  are  obtained.  These  results  are 
useful  in  robustness  studies  of  tolerance  limits, 
estimates  and  statistical  tests  derived  for  specified 
distributions  such  as  the  normal  distribution.  Some 
examples  are  given. 


VT 
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TOLERANCE  LIMITS  AND  INEQUALITIES  FOR  s-ORDERED  DISTRIBUTIONS 

by 


Michael  J.  Lawrence 


I.  INTRODUCTION 

Motivated  by  classical  but  somewhat  unsatisfactory  measures  of  skewness 
and  kurtosis,  van  Zwet  (1964)  introduced  two  partial  orderings  on  the  space  of 
distribution  functions.  Barlow  and  Proschan  (1966a)  investigated  the 
properties  of  linear  combinations  of  order  statistics  from  distributions 
ordered  in  the  sense  of  van  Zwet  but  corresponding  to  positive  random  variables. 
Our  objective  is  to  extend  the  results  of  Barlow  and  Proschan  (1966a, 1966b)  to 
distributions  s-ordered  in  the  sense  of  van  Zwet  but  not  restricted  to  the 
positive  axis.  If  we  were  to  confine  attention  exclusively  to  symmetric 
distributions,  this  would  be  a  relatively  straightforward  task.  However,  we 
extend  the  s-ordering  definition  of  van  Zwet  to  include  a  wider  class  of 
possibly  skewed  distributions.  We  obtain  tolerance  limits  which  are 
conservative  for  a  wide  class  of  distributions.  This  class  of  distributions 
is  especially  reasonable  for  measurement  type  data.  In  addition,  we  obtain 
many  results  concerning  linear  combinations  of  order  statistics  from  s-ordered 
distributions.  These  results  should  be  useful  in  robustness  studies  of 
estimates  and  statistical  tests  derived  for  specified  distributions  such  as 
the  normal  distribution. 

The  basis  of  van  Zwet's  ordering  between  distribution  functions,  and 

hence  between  random  variables,  is  that  one  random  variable  can  be  expressed 

as  a  convex  or  concave-convex  transformation  of  another  random  variable.  We 

adopt  van  Zwet's  notation  for  c-ordering:  F  <  G  if  and  only  if  G  ^F  is 

convex  on  the  support  of  F  ;  and  a  more  general  definition  of  s-ordering: 

F  <  G  if  and  only  if  F(m)  ■  G(m)  ■  \  and  G  ^F  is  concave-convex  about  m  , 
s 


on  the  support  of  F  .  If  F  and  G  are  symmetric,  then  our  definition  of 
s-ordering  coincides  with  that  of  van  Zwet.  For  convenience,  we  shall  assume 
that  the  median  of  F  and  G  is  the  origin. 

One  expects  that  If  F  <  G  then  G  Is  more  skewed  to  the  left  than  F  . 
If  skewness  is  measured  by  the  standardized  odd  central  moments,  then  we  would 
expect  (assuming  for  convenience  that  EX  ■  EY  ■  0  ,  and  that  X  (Y)  has 
distribution  F  (G)  ) 


(1.1) 


EX2^1  <  EY2^1 

j  k+4  —  .  k+^s 

(EX  ;  (EY  ; 


Van  Zwet  proves  this  result  and  also  that  if  i  and  n  tend  to  infinity  and 

lim  ■^■r,0<r<l,  then  F  <  G  if  (1.1)  holds  asymptotically  for  some 

fixed  k  ,  all  0  <  r  <  1  ,  and  order  statistics  X.  and  Y .  .  If  F 

i,n  i,n 

and  G  are  symmetric,  then  F  <  G  would  seem  to  imply  that  G  has  heavier 

s 

tails  or  more  "kurtosis"  than  F  .  Van  Zwet  proves  this  when  the  even 

standardized  central  moments  are  taken  as  the  measure  of  "kurtosis". 

Barlow  and  Proschan  (1966a,b)  have  derived  tolerance  limits  for  the 

distributions  which  are  c-ordered  with  respect  to  certain  distribution 

functions,  as  well  as  developing  many  interesting  inequalities  for  linear 

combinations  of  the  order  statistics  from  c-ordered  distributions.  They  have 

particularly  exploited  the  properties  of  distribution  functions  which  are 

convex  with  respect  to  the  exponential  distribution,  and  have  also  introduced 

and  developed  for  positive  random  variables,  the  properties  of  the  weaker 

G~^F(x) 

star-shaped  ordering;  i.e.,  - J  is  increasing  in  x  for  x  on  the 

support  of  F  and  F(0)  ■  0  . 


It  is  the  object  of  this  thesis  to  develop  some  of  the  statistical 
properties  of  distribution  functions  related  by  both  s-ordering  and  star- 
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shaped  ordering.  We  will  be  particularly  interested  in  a  natural  class  of 
distribution  functions — viz.,  those  which  are  s-ordered  with  respect  to  the 
Laplace  distribution  (also  called  the  bilateral  exponential  distribution). 

Examples  of  such  s-ordered  distributions  are: 

U-shaped  <  uniform  <  normal  <  logistic  <  Laplace  <  Cauchy 

8  8  8  S  S 

(cf.  Van  Zwet  (1964),  pp.  70-71,  72-73.) 

There  are  numerous  reasons  why  we  are  interested  in  studying  the  properties 
of  distributions  related  by  s-ordering.  The  most  apparent  one  is  that  very 
often  we  do  not  know  the  exact  distribution  but  because  of  physical  consider¬ 
ations  we  can  make  certain  deductions  about  the  properties  of  the  distribution, 
and  hence  do  not  want  the  disadvantage  of  a  distribution-free  approach.  For 
instance  if  we  measure  the  length  of  an  object,  it  is  plausible  that  the 
probability  of  obtaining  an  error  in  the  range  (|x|,|x|  +  5  x)  ,  given  that  the 
error  is  at  least  |x]  is  increasing  in  |x|  .  The  class  of  distributions 
having  this  property  are  s-ordered  with  respect  to  the  Laplace  distribution. 
Another  example  would  be  the  commonly  occurring  situation  where  the  normal 
distribution  is  assumed,  but  we  suspect  that  this  is  not  true  and  that  in  fact 
the  tails  are  heavier  or  lighter  than  the  tails  of  the  normal.  We  then  wish 
to  know  if  the  normal  assumption  is  conservative  or  not.  In  short,  we  want  to 
know  how  robust  the  normal  distribution  is  against  s-ordered  alternatives. 

In  Chapter  II,  we  develop  some  new  inequalities  for  concave-convex 
functions  <£  ,  when  $  satisfies  a  skewness  condition.  These  Inequalities 
are  used  in  Chapter  III  to  construct  tolerance  limits  for  distributions  s- 
ordered  with  respect  to  the  Laplace  distribution  (these  distributions  we  call 
SIFR).  We  further  investigate  symmetric  SIFR  distributions  by  obtaining  sharp 
bounds  on  F  symmetric  and  SIFR  when  given  only  the  mean  and  the  variance  of 
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F  .  This  is  used  to  construct  a  confidence  bound  on  the  variance  of  F  . 

Also  in  Chapter  III,  the  robustness  of  tolerance  limits  from  the  normal 
distribution  against  s-ordered  alternatives  is  investigated.  Results  are 
obtained  for  small  sample  sizes. 

In  Chapter  IV,  we  are  interested  in  inequalities  for  the  order  statistics 

and  their  expectation  when  related  by  the  weaker  r-ordering.  We  say  that 

F  <  G  if  F(0)  ■  G(0)  -  and  — — is  increasing  (decreasing)  in  x 

positive  (negative).  Bounds  on  the  expectation  of  the  order  statistic 

from  F  are  given  in  terms  of  the  expectation  of  the  order  statistics  from 

G  when  F  <  G  and  F  and  G  are  symmetric.  Also,  an  Inequality  relating 

a  linear  combination  of  the  expectations  of  the  order  statistics  from  F  and 

G  is  given  for  F  <  G  ,  G  symmetric  about  the  origin  and  the  direction  of 

s 

the  skew  of  F  known. 

If  F  <  G  and  F  and  G  are  symmetric,  then  we  prove  that  not  only  are 
the  standardized  even  central  moments  of  G  greater  than  that  of  F  but  so 
are  the  usual  sample  estimates  of  the  standardized  even  central  moments  in  a 
stochastic  sense. 

Preliminaries 

We  adopt  the  following  definitions: 

(i)  F  <  G  if  and  only  if  G  *"F  is  convex  on  the  support  of  F  . 

(ii)  F  <  G  if  and  only  if  F(0)  ■  G(0)  -  4  and  G  ^F  is  concave-convex, 

8 

about  the  origin,  on  the  support  of  F  . 

(iii)  F  <  G  if  and  only  if  F(0)  -  G(0)  -  4  and  ^  *s  *-ncreas*n8 

(decreasing)  for  x  positive  (negative)  on  the  support  of  F  . 

(iv)  F  is  SIFR(SDFR)  if  and  only  if  F  <  (>)  G  ,  where  G  is  the  Laplace 

8  S 

-Ixl 

distribution;  i.e.,  G' (x)  “  e  1  1  ,-°°<x<°°. 

(v)  X  Y  if  and  only  if  P(X  <_  a)  >_  P(Y  <_  a)  ,  <  a  <  «  . 


V 


We  assume  throughout  that  F  and  G  are  continuous  with  median  at  the 

origin.  Mote  that  F  <  G  implies  F  <  G  . 

s  r 

Define  F^  to  be  the  distribution  of  the  i-th  order  statistic  from  F 


and  F  to  be  the  distribution  of 
a 


Since  G^F(i) 


--I* 

G  F  we -see 


(a)  F  <  G  implies  F(i)  <  G{i) 

(b)  F  <  G  implies  F,.v  <  G/JX 

r  Y  (i)  r  (i) 

If  furthermore,  F  and  G  are  symmetric  about  the  origin,  then  for  o  >_  1 


(c)  F  <  G  implies  F  <  G 
s  a  c  a 


(d)  F  <  G  implies  F  <  G 
r  a  r  a 


on  the  positive  axis. 

If  G  is  symmetric  about  the  origin  and  Jb  ■  (G(0x)  |  0  >  0}  then  a 

sufficient  statistic  for  £  based  on  a  complete  random  sample 

*Y1,Y2 . Yn)  is  8*ven  by  ( | Yx | ,  |Y2 1 ,  ...,  |yJ)  .  Suppose  we  are 

interested  in  studying  the  robustness  of  statistics  derived  under  the 

assumption  that  the  observations  are  distributed  according  to  G  when  in  fact 

they  are  distributed  according  to  F  where  F  <  G  .  Since  by  (c) ,  F  <  G 

8  8 

implies  F^  <  G^  ,  the  results  of  Barlow  and  Proschan  (1966a)  apply  to  linear 

combinations  of  the  sufficient  statistics  for  G  . 

Throughout  we  let  X,  <  X„  <  . . .  <  X  (Y,  <  Y„  <  . . .  <  Y  )  be  an  ordered 

l  —  i.  —  —  n  i.  ~  L  —  —  n 

sample  from  F  (G)  ,  and  we  observe  that  Y  G  *F(X)  where  denotes 

stochastic  equality. 


i 
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II.  INEQUALITIES  FOR  CONCAVE-CONVEX  FUNCTIONS 

Some  inequalities  on  concave-convex  functions,  interesting  in  their  own 
right  and  also  required  later  on,  will  now  be  developed.  We  will  need  a  slight 
extension  of  a  theorem  by  Hardy,  Lit tlewood  and  P<51ya  (1929).  In  this 
connection,  see  also  Barlow,  Marshall  and  Proschan  (1967)  and  Karlin  and 
Novikoff  (1963)  p.  1252. 

We  say  that  $  is  concave-convex  about  the  origin  and  defined  on  [-a,b] 
if  $  is  concave  on  [-a,0]  and  convex  on  [0,b]  . 


Theorem  2.1 

Let  y  be  a  signed  measure  on  [-a,b]  ,  0  <  a  ,  b  <  °°  ,  then 


(2.1) 


b 

/ 

-a 


4>(x)dy(x)  >  0 


for  all  <p  concave-convex  about  the  origin,  continuous  at  the  origin,  and 
defined  on  [-a,b]  if  and  only  if 


(2.2) 


b  b 

xdy(x)  ■  0  ,  J" du(x)  ■  0 
-a  -a 


and 


-z '  b 

(2.3)  J " (x  +  z')du(x)  +  J  (x  -  z)dy(x)  0 

-a  z 


for  all 


z  e  (0,b] 
-z '  c  (-a,0] 
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Proof 

Suppose  first  that  inequality  (2.1)  is  satisfied  for  all  <p  concave- 
convex  about  the  origin,  continuous  at  the  origin  and  defined  on  [-a,b]  . 
Now  (2.2)  follows  if  <j>(x)  ■  (+)  1  or  cp (x)  *  (+)x  ,  and  (2.3)  follows  if 
<j>(x)  is  a  "double  angle  function,"  i.e.,  a  function  of  the  form 


Next,  suppose  that  (2.2)  and  (2.3)  hold.  Since  <p  is  concave-convex 
about  the  origin  and  continuous  there,  there  exists  a  ,  8  such  that 

4>(x)  -  ax  -  8  0  for  x  >_(<)  0 

Hence  we  may  assume  that 

<Hx)  2.(<)  0  for  x  >_(<)  0  . 

Consider  now  a  sequence  of  functions  $n(x)  such  that  each  $  is  the 
sum  of  a  finite  number  of  positive  multiples  of  "double  angle  functions"  and 
d>n (x)  is  increasing  in  n  for  x  ^  0  and  decreasing  in  n  for  x  <_  0  . 

Since  <p  is  concave-convex  we  can  construct  a  sequence  of  functions 

<t>n(x)  >_(<)  0  for  x  >_(<)  0  such  that  $n(x)  converges  upward  (downward)  to 

$  for  x  ^  0  (x  <  n)  .  By  construction 

b 

j m $n(x)dy(x)  >_  0  , 

a 

for  each  and  (2.1)  follows  from  the  Lebesgue  monotone  convergence  theorem.  | | 


r 
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Remark 

If  $  in  Theorem  2.1  is  not  only  concave-convex  but  also  satisfies 
4>(x)  <_  -<K-x)  ,  then  we  see  that  the  inequality  (2.3)  need  only  hold  for  all 
z  ,  z'  satisfying  z  z'  . 

We  will  now  prove  a  lemma  which  although  fairly  elementary  is  nevertheless 
quite  useful. 


Lemma  2. 2 


n 


l 

i-1 


Vi 


+  l  b  6  >  x 

j-1  J  J 


for  all 


1  ilf,  >  Y,  >  ....  >  y  >0 
—  1—  2—  —  n  — 

1  6  >  6  >....>  6  >  0 
1  —  2  —  —  m  — 

if  and  only  if 


1  k 

1  ai  +  l  bi  1  x  for  0  <_  j  <_  n 

1  1  and  0  k  £  m 

0 

and  we  define  £  -  0  . 

1 

Proof 

Let 


u 

o 


n 


u . 
J 


'j  Tj+1 


u 

n 


Y 


n 


v  * 


o 


v  *  <5 

m  m 


and 


9 


A 


j 


The  lemma  can  now  be  reformulated  as 


(2.4) 


m 


l  Vi + 1 


Vk  i x 


for  all 


0  <  u  <  1  ,  0<j<n 

-  j  -  ~  J  - 

(2.5) 

0£vk_<l  ,  0  <  k  <  m 


n 

and  l  u 
0  J 

m 

and 


1 


1 


if  and  only  if 


(2.6) 


+  B  >  x 
k  — 


for  0  <_  j  _<  n 
0  <  k  < 


The  proof  is  straightforward.  If  (2.4)  is  true  for  all  u^  , 
(2.5),  then  clearly  the  inequality  (2.6)  must  be  satisfied, 
then 


l 


k-0 


x  , 


which  completes  the  proof.  | | 

We  would  like  to  determine  conditions  on  a^,a2»  . ..,  a^ 


n 

<_  [  a.$(x.) 
1 


holds  for  all  <_  X2  ^ 


<  x  and  all  $  concave-convex 
—  n 


satisfying 
If  (2.6)  is  true, 


such  that 


.  However,  it  is 
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possible  to  construct  an  example  to  show  that  no  such  Inequality  holds  for  all 


$  concave-convex,  all  ordered  x's  and  a^  0  ,  1  <_  i  <  n  .  If  we  assume 
that  Is  concave-convex  about  the  origin  and  $(0)  ■  0  ,  then  a  simple 
additional  condition  on  <J>  which  admits  a  solution  to  our  problem  is 
<Kx)  _<  -<K-x)  for  x  ^  0  . 

In  proving  the  next  theorem  it  is  necessary  only  to  consider  <p  such 
that  $(0)  *  0  ,  and  <j>  is  concave-convex  about  the  origin.  This  can  be  seen 


by  making  the  transformation  $  (x)  ■  <Kx-*-c)  -  <Kc)  ,  if  d>  Is  concave-convex 


about  x  ■  c  .  To  simplify  notation  we  will  define  A 

**■!■*• 


j  ■ !  ai 


Theorem  2.3 


If  <J>  la  concave-convex  about  the  origin  and  defined  on  (-°°,b)  , 

<Hx)  <  — <K-x)  for  x  _>  0  ,  4> (0)  *  0  and  <t>  is  continuous  at  the  origin  then 


(2.7) 


Ai  Vl)  1 1  ‘v*<Xli 


for  all  Xj  <  *2  £  •••  _<  x^  _<  0  <_  Xj^  <_  . . .  <_  xn  <_  b  if  and  only  if  (2.8) 


or  (2.9)  is  satisfied,  where 


(2.8) 


0  1  -Ai  +  Aj  1  1 
At  1  0 
^  L  0 


for  l<^i<^k<j<_n  , 


A  -  A  -  0 
q  p 

(2.9)  l<Ar,<  Ar,+1  1  .. 

A  <■  —  1 

Aj  -  1 


.  <_  A,  >  for  0<p<r'_f_ifk<j<_r<q<n  , 


and  r'  (1  <_  r'  <_  k  +  1)  ,  r(k  <_  r  <_  n)  are  fixed. 
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Proof 

It  will  be  convenient  to  assume  ■  0  with  a^  ■  0  .  Clearly  we  can 
always  do  this  with  no  loss  of  generality. 

We  use  Theorem  2.1  and  adopt  the  measure 

n 

-  1  ,  x  -  l  ax 

1 

aA  ,  x  -  xt 

n 

1  -  I  a.  ,  x  -  0 
1 

0  ,  elsewhere 

Clearly  then 


oo 


du  (x) 


00 


xdu(x) 


0 


0  .  Similarly 


12 


and 


0  i  \_i  i  1  <  xt  • 


If  no  x  exists,  add  an  extra  term  x  <  -z'  and  a 

8  0  0 


add  x  t  >  z  if  x  <  z  and  let  a  ,  »  0  . 

n+l  n  —  n+1 


Case  (1) 


n+1 


Assume  that  -z'  \  a.x  <_  z  .  Now  the  left  hand  side  of  (2.10)  equals 

0  1 


(2.11) 


8-1 

-  I  A1(x1  -  x1+1)  -  4a(xa  +  X')  -  At(xt  -  2)  - 
n+1 

-  I  A  (x  -  x  )  • 

J.+1  1  1 


Since  x^^  -  <_  0  ,  0  and  since  these  differences  can  be 
arbitrarily  small  for  suitable  choice  of  the  x  values,  we  see  that  the 
necessary  and  sufficient  conditions  for  (2.11)  to  be  nonpositive  are 


<  0  H  i  <k 

(2.12) 

_>  0  k  <  j  <_  n 


If  (2.12)  holds  then  )  a.x.  >  0  . 

*•  i  i  — 


Case  (2) 

n 

Assume  z  <  J  a  x.  ,  and  recall  that  z  ^  z'  >_  0  .  The  left  hand  side  of 
1  1 

(2.10)  can  be  written  as 


13 


k-1  n+1 

l  Ai(xi  -  xi+1)  +  I  A1(x1 


k+1 


xi-i) 


-  z 


8-1 

'  A8<XB  +  Z>)  -  I  Ai(Xi 


Xi+1) 


n+l 

-  A  (x  -  z)  -  l  A  (x 
£  £  £+1  1  1 


Xl-1) 


(2.13) 


71  (_AS  +  W  +  •  * 


x£,_l  _  _  X£-2 

+  — r  (-A* +  Ai-i> +  z  (-At-i +  At-2)  + 


.  z  — 

Now  from  Lemma  2.2  and  identifying  x  *  -  A  +  A  -  1  ,  we  see  that  the 

Z  S  A* 

necessary  and  sufficient  conditions  for  (2.13)  to  be  nonpositiva  are 


—  A  +  A  -  1  <  0 
z  s  £  — 


(2.14) 


for  8<h<^k<g<_£ 


+  A  -1-Ah<_0 


Since  (2.14)  must  hold  for  all  z  ,  z'  such  that 


z'  <  z  < 


l  Vj  . 


then  from  (2.12)  and  (2.14), 


-  —  A  +  A  -  1  <  -A  +  A  -  1  , 

z  s  £  —  s  £ 


M) 


+  A  -  1  -  A,  <  A  -  1  -  A. 
e  n  -  R  n 


and  we  have  that  (2.14)  holds  if  and  only  if  -A,  +  A  <1  for  0  <  h  <  k  , 

h  g  -  —  — 

k  <  g  <  n  . 


Thus,  we  have  established  the  necessity  and  sufficiency  of  conditions 
n 

(2.8)  for  the  case  £  a,x,  >_  -z'  . 

J-l  1 


Cage  (3) 

Assume  that  T  a.x.  <  -z'  .  Note  that  if  -z'  <  x,  and  z  >  x  then 
Mi-  —  1  —  n 

the  left  hand  side  of  (2.10)  is  nonpositive.  The  left  hand  side  of  (2.10) 
can  be  written  as 


(2.15) 


:  [/a.  -  —  A  +  ^) 
y  SL  z  s  z  f 


s+1 

z 


(-A, 


+  W 


+ 


(-Afl 


+  kl-l> 


By  Lemma  2.2,  the  necessary  and  sufficient  conditions  for  (2.15)  to  be 
nonpositive  are  that 


A 


l 


for  s 


h  £  k  <  g  <  l 


—  (1-A  )  +  A  -  A,  +  A  <0) 
z  s  s  h  g  —  / 

By  noting  that  if  -z'  <  x.  then  A  ■  A„  *  0  ,  and  if  z  >  x  then 

A.  -  A  . ,  -  0  ,  and  that  the  above  conditions  must  hold  for  all  z  >  z'  , 
l  n+1  — 

we  see  that  the  necessary  and  sufficient  conditions  for  (2.7)  to  be  true  in 
the  case  that  £  aixi  ~  ~z'  *  are  (2.9).  Now  from  (2.9)  we  see  that 
\  a^x^  x^  and  hence  if  \  aixi  il  -z'  then  the  left  hand  side  of  (2.10) 
is  zero.  |  | 


Remarks 

1.  We  see  that  the  only  concave-convex  functions  which  admit  a  solution 

of  (2.7)  with  all  weights  a^^,  nonzero  are  those  which  can  be 

z  * 

generated  from  a  double  angle  function  with  —  bounded  above  and  below. 

2.  If  in  Theorem  2.3,  we  reverse  the  skewness  condition  on  $  and  set 


$(x)  2.  -<K-x)  for  x  >_  0  ,  then  we  can  see  from  the  proof  of  the  theorem 


that  (2.8)  is  replaced  by 


0  -  *j  -  1 


Ai  ■  0 


for  l<^i<^k<j<^n  , 


and  for  (2.9)  we  need  r'  ■  k+1  . 

By  considering  d>  (x)  ■  -4>(-x)  we  get  the  following  corollary. 


Corollary  2,4 

If  concave-convex  about  the  origin  and  defined  on  (-a,00)  ,  $(0)  ■  0 
4)  continuous  at  the  origin  and  4>(x)  ^-$(-x)  for  all  x  >_  0  ,  then 


*(|  Vi)  -  \ 


for  all  -a  x1  ^  x2  _<  . . .  ^  xk  _<  0  xfe+1  <  ...  <  xfl  if  and  only  if  (2.17) 


or  (2.18)  is  satisfied. 


(2.17) 


Ai  >0 


Aj  -  ° 


0  1  Ai  -  Aj  I  1 


for  1  <  i  <  k  <  J  <n  , 


A  ■  A  *  0 

q  p 


(2.18)  1  <  Af  <  Af+1  1  ...  1  Afe+1  >  for  0<p<r'<i<_k<J_<r<q<n 


A,  1-1 


and  r'  (1  <_  r'  <_  k+1)  ,  r(k  <^r  <  n)  are  fixed. 


r 
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By  noting  that  if  <J>  la  concave-convex  then  -<p  is  convex-concave,  we 
can  obtain  similar  inequalities  for  convex-concave  functions.  As  they  follow 
trivially,  we  shall  not  include  them  here. 

In  the  special  case  that  <p  is  concave-convex  and  antisymmetric  (i.e., 
$(x)  ■  -<K-x))  i  we  can  set  z  ■  z*  in  the  proof  of  Theorem  2.3  and  obtain 
the  necessary  and  sufficient  conditions  on  a^^,  •••»  an  such  that  (2.7) 
is  true  for  all  <J>  concave-convex  and  antisymmetric.  However  for  this  special 
case  the  result  also  follows  directly  from  Lemmata  4.1  and  4.3  of  barlow  and 
Proschan  (1966) . 

Theorem  2.5 

If  is  concave-convex  and  antisymmetric  about  the  origin  and  defined 
on  (-00,00)  then 


*(!  vi)  i  |  ai*(xi> 


for  all  x^  <_  . . .  <_  x^  0  .  _<  i  If  and  only  if  (2.8)  or  (2.19)  is 

satisfied. 


(2.8) 


0  1  -ki  +  Aj  1  1 


A.  <  0 


i  - 

>1 

C— u 

1  v 

O 

Ai  -  1 

) 

2  for 

for  1<  i  <  k  <  j  <_  n 


A  *  A 


0 


■>*>»<&  '  i"™'  •• 


Proof 


The  proof  follows  from  Barlow  and  Proschan  (1966)  with  the  observation 


ns  n 

l  a14»(x.)  -  l  -  a  4>(-x  )  +  £  a.*(x.)  , 

1  1  1  1  k+1  1  1 


and  hence  the  right  hand  side  only  involves  $  convex  on  [0,®)  . 


Theorem  2.6 

If  $  is  concave-convex  and  antisymmetric  about  the  origin  and  defined 


on  then 


for  all  ^  *2  —  •  •  *  -1  xk  —  ®  —  •  •  •  —  xn  and  only  if  either  (2.20)  or 
(2.21)  is  satisfied. 


(2.20) 


Ai  I  0 


Aj  I  0 


0  1  Ai  -  Aj  I  1 


for  1  1  1  1  k  <  j  <  n 


Ail-1 


(2.21)  A  >  1 


A  ■  A  ■  0 
p  q 


for  0  <  p  <  r'  liik<  jir  <  q  <n  . 
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III.  TOLERANCE  LIMITS  AND  CONFIDENCE  BOUNDS 

It  will  be  convenient  to  let  X  (Y)  have  distribution  F  (G)  .  We 
assume  that  G  is  symmetric  about  the  origin  and  strictly  increasing  on  its 
support  and  that  both  F  and  G  are  continuous.  Let  X,  <  X0  <  . . .  <  X 
(Y^  <  Yj  ^  . . .  <_  Y  )  denote  an  ordered  sample  from  F  (G)  .  We  say  that  a 
random  variable  X  is  stochastically  greater  than  a  random  variable  Y  , 
denoted  by  X  Y  ,  if  P(X  _>  x)  >  P(Y  ^  x)  for  all  x  . 

We  construct  tolerance  limits  by  using  the  relevant  inequalities  of 
Chapter  II  to  give  stochastic  comparisons  between  the  order  statistics  from 
F  and  G  when  F  <  G  .  Consider  for  example  Theorem  2.3.  The  weights 

S 

a^.a^i  . . . ,  an  in  inequality  (2.7)  which  must  be  selected  to  satisfy 
conditions  (2.8)  or  (2.9),  are  dependent  on  the  value  of  the  index  k  defined 
by  X^  <_  .  . .  <_  X^  <_  0  X^+1  <^  . . .  <_  X^  ,  which  is  therefore  a  random  variable. 

Let  X  «  (X^,X2,  . Xn)  denote  a  random  vector  of  ordered  observations. 

Thus  when  making  a  stochastic  comparison  using  inequality  (2.7),  the  weights 

must  be  chosen  as  a  function  of  X  ,  i.e.,  (a^(X) . an(X))  =  [a(X)  ]  . 

We  require  that  these  weights  satisfy  conditions  (2.8)  or  (2.9)  for  every 
possible  outcome,  X  ;  and  then  say  that  [a]  satisfies  (2.8)  or  (2.9). 

If  we  let  <J>(x)  ■  G  ^F(x)  in  Theorem  2.3,  then  the  condition 
G  ^F(x)  <_  -G  ^F(-x)  ,  x  ^  0  ,  will  be  satisfied,  for  example,  if 
G(0)  *  F(0)  ■  ,  G  is  symmetric  about  the  origin  and  1  -  F(x)  >_  F(-x)  for 

x  >_  0  . 

Theorem  3.1 

If  F  <  G  ,  G  symmetric  about  the  origin,  F(0)  =*  G(0)  m  h  , 
s 

1  -  F(x)  >  F(-x)  for  x  _>  0  and  [^]  satisfies  (2.8)  or  (2.9),  then 


G(!ai®Yi)  ?  f(\  si®xi)  • 
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Proof 

From  Theorem  2.3  we  have  that  if  [aj  satisfies  (2.8)  or  (2.9),  then 


The  stochastic  equality  follows  from  the  fact  that  G  preserves 
order  with  respect  to  the  origin,  and  G  ^F(X^) ,  G  ^F(X2) ,  ...,  G  ^F(Xn)  are 
jointly  distributed  as  the  order  statistics  from  G  .  | | 

Corollary  3.2 

If  F  <  G  ,  G  symmetric  about  the  origin  F(0)  -  G(0)  -  , 

s 

1  -  F(x)  <_  F(-x)  for  x  ^0  and  [a]  satisfies  (2.17)  or  (2.18),  then 

G(|*i®Yi)  1*  . 

If  we  let  F  be  symmetric  about  the  origin,  we  can  similarly  prove  the 
following: 

Theorem  3.3 

If  F  <  G  ,  F  and  G  symmetric  about  the  origin  and  [aj  satisfies 
s 

either  (2.8)  or  (2,19),  then 


G(|ai®Yi)  -  F(\ ai®xi 
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Corollary  3.4 


If  F  <  G  ,  F  and 
either  (2.20)  or  (2.21) , 


G  symmetric  about  the  origin,  and 
then 


[.a]  satisfies 


G 
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3.1  Tolerance  Limits  for  SIFR  and  SDFR  Distributions 

Barlow  and  Proschan  (1966b)  have  considered  tolerance  limits  for  the 
class  of  distributions  which  are  c-ordered  with  respect  to  the  exponential. 
These  distributions  possess  an  increasing  failure  rate,  and  so  arise  naturally 
when  wear-out  is  present.  The  two  sided  analogue  of  the  increasing  failure 

rate  class  is  the  class  of  distributions  which  are  s-ordered  with  respect  to 

-Ixl 

the  Laplace  distribution  G  ,  i.e.,  G'  (x)  ■  %  e  1  1  .  We  will  call  the 

distribution  F  SIFR(SDFR)  when  F  <  G  ,  (F  >  G)  and  G  is  the  Laplace.  We 

s  s 

note  that  this  implies  that  F(0)  *  G(0)  ■  ^  .  The  distribution  function  F 

with  density  f  ,  F(0)  *  is  SIFR(SDFR)  if  ^  *S  *ncreasin8 

f  (x) 

(decreasing)  in  x  for  x  ^  0  and  y(x)  decreasing  (increasing)  in  x 

for  x  <_  0  .  Note  that  F  need  not  be  symmetric  and  that  we  chose  the  median 

at  the  origin  only  for  convenience.  Recall  that  the  normal  distribution  is 

SIFR  and  the  Cauchy  distribution  is  SDFR. 

The  SIFR  class  can  arise  naturally  when  we  are  considering  problems  such 

as  the  distribution  of  the  error  of  some  measurement,  for  here  we  would  expect 

that  the  probability  of  an  error  in  the  range  (|x|,|x|  +  6x)  ,  given  the 

error  is  greater  than  |x|  ,  would  be  increasing  in  |x|  .  We  will  develop 

tolerance  limits  for  F  SIFR  and  F  SDFR  in  the  one  sided  case  when  F  is 

skewed  and  in  the  two  sided  case  when  F  is  symmetric. 

If  we  have  a  complete  sample  and  F  is  symmetric  about  the  origin,  then 

it  is  an  easy  matter  to  construct  conservative  two  sided  tolerance  limits  for 

\  —|x|/0 

F  .  We  need  only  assume  F  <  G  where  G'(x)  ■  rr  e  .  Then 

r  zt ) 

(IyL  <  IyL  <  ...  <  I Y I  )  is  a  sufficient  statistic  for  0  and  tolerance 

limits  for  G  may  be  constructed  from  these.  Since  F^  ,  the  distribution 

of  | X |  is  star  shaped  with  respect  to  G  ,  the  distribution  of  |y| 

3  G_1Ffl(x) 

(which  is  the  exponential  distribution)  i.e.,  -  is  increasing  in 
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x  ,  we  can  apply  Theorem  3.1  of  Barlow  and  Proschan  (1966b).  Let 


and 


er>nQ0  -  l  (n  -  1+1) r'1  ( (X|  ±  -  |X|  i_1) 


B  «  ~2r  log(l-g) 

ct.q.r  2.„  v 

Xa(2r) 


and 


Then 


** 

C  *=  max  (B  .  r(n  -  r+1)  ) 

a,q,r  a,q,r’ 


**  * 

C  0  (X) 

a,q,r  r,n  — 


** 

-c  e  (x) 

a,q,r  r,nv- 


dF(x)  >_  q  \  >_  1-a 


However,  If  we  have  a  censored  sample  and/or  F  is  not  symmetric  this 
inequality  is  not  valid. 

We  will  need  to  develop  some  properties  of  the  Laplace  distribution. 

Given  an  ordered  sample  X,  <  X„  <  . . .  <  X.  <  0  <  X,  .<...<  X  ,  we  define 

1-2-  —  k  —  —  k+1  —  -  n 


the  statistic  0  .  (X)  by 

r  ,  r ,  n  —  J 


,r,n®  '  l,  -  1<Xi  -  W  +  j2  <"  -  1+1)<Xi  -  Xi-1> 


If  X,  >  0  let  k  -  1  and  if  X  <  0  let  k  =  n-1  . 
1  n 
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We  can  easily  verify  that  if  <_  X2  <_  . .  .  £  X^  £  0  £  X^+1  <_  ...  Xr 
are  distributed  as  the  n  order  statistics  from  distribution  G  where 
G'  (x)  *  —  e  I  xl^e  t  Q  >  o  (the  Laplace  distribution  with  scale  parameter  0  ), 
then  the  maximum  likelihood  estimate  of  the  scale  parameter  0  ,  given  censor¬ 
ship  at  r'  and  r  where  r'  <  k  ,  r  >  k  is 


-Xr,(r'-1)  +  I,  |X±|  +  Xr(n-r) 


r-r'+l 


r "  it 

Hence  we  see  that  if  r  and  r'  are  such  that  lim  —  and  lim  —  exist 

0  ,  (X)  n  n 

r  *  r  n  — 

finitely  then  — —  is  asymptotically  the  maximum  likelihood  estimate 
for  the  scale  parameter  when  X  has  the  Laplace  distribution  with  scale 
parameter  0  . 


Lemma  3.5 


If  Y  ~  G  where  G'  (y)  *  h  e 


’  fvl  =  U  p  ^  y 


then 
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r,  r  n(D  '  x  (2(r-r#-l))  . 

L  |  L  f  11 


(N.B.  denotes  "is  distributed  as".) 


Proof 


If  Y  -  G  ,  the  Laplace  distributed  function,  then 


P(Y  <  x  +  y  |  Y  >  y  >  0)  «  1  -  e 


-x 


which  is  independent  of  y  .  Now  if  Y^  is  the  hC^  order  statistic  from 
G  we  have 


p(VrYh-xi  Vy-0)  "  p(Yh+i-x+yiYh“y-0)  =  e 


-(n-h)x 
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and  since  this  is  independent  of  y 


P<Vl  '  Vh  >  x  |  Yh  >  0) 


- (n-h)x 
e 


or  the  random  variable 


<”-h>(Vl  -  V 


given  that  Y^  _>  0  is  distributed  as  the  unit  exponential.  Now  from 
symmetry  we  have 


p  (Y  -Y  u^>x|Y  _.<0) 
n-h  n-h+1—  1  n-h+1— 


p(Yhfi-VxIV°) 


- (n-h)x 


and  hence  the  conditional  distribution  of  i(Yi+1  -  Y^  given  Y^+1  <_ 
exponential  with  mean  unity,  which  completes  the  proof.  |  | 

Let  H  (r',r)  -  ,  and 

“*q  X„2(r-r'-l) 

a 

if  x^2(r-r'-l)  ^  -2m  log(l-q) 
if  x^2(r-r'-l)  >  -2m  log(l-q) 


H  (r * ,r) 
a  ,q 


C  n(r’’r)  " 
a ,  q 


l 

m 


where  m  »  min(r' ,n-r+l)  . 

We  first  give  the  tolerance  limit  on  the  tail  probability,  1  -  F 
values  to  the  left  of  the  median. 


Theorem  3.6 

If  F  is  SIFR,  F (0)  -4,1-  F (x)  >_  F(-x)  for  x  >_  0  ,  and 
0  q  <_  4  then 


Pf{  1  -  FK,I-2q(r’-r>V,r,n<«]  1-0  • 


0  is 


,  for 


25 


Proof 


Let  G  be  the  Laplace  distribution,  i.e.,  G'(x)  ■  *s  e  '  ,  then 

from  Theorem  3.1, 


f[-H  .  _  (r',r)6  ,  (X)l  g|"-H  ,  _  (r',r)6  ,  (Y)l 

L  a » l~2q  r  , r , n  —  J  -  [  a,l-2q"  *  r',r,nv-'J 


for  mH  ,  _  (r',r)  >  1  .  It  follows  in  this  case  that  since 
a,l-2q  - 


MGK.l-2,(r'>r)5r'.r,„<X>] 


we  have 


PF{FK,l-2<j(r,'r)5r',r,„®J  -l)  i1- 


Now  if  mH  1  0  (r',r)  <  1  ,  we  note  that 
a i i-^q 


F  [  -  -  0  ,  (X)\  <l  G  (-  -  6  ,  ( Y ) J  <  G  [  -H  .  .  (r',r)6  ,  (Y)) 

\  m  r  ,r,n  —  J  -  \  m  r',r,nv-y  -  ^  a,l-2q  ’  r',r,n-y 


which  proves  the  theorem.  | | 

Clearly  Theorem  3.6  implies  a  tolerance  limit  on  the  distribution  F  for 

values  to  the  right  of  the  median.  However,  we  shall  not  present  this  result. 

We  note  that  the  tolerance  limit  in  Theorem  3.6  achieves  1-a 

confidence  for  the  Laplace  distribution,  (i.e.,  is  sharp)  only  when 
2 

Xa [ 2 (r— r  — 1 ) ]  <_  -2m  log2q  .  The  following  corollary  gives  bounds  on  q  and 
a  such  that  this  condition  is  satisfied. 


Corollary  3. 7 


If  F  is  SIFR,  F (0)  =*5,1-  F(a)  >_  F(-x)  for  x  >_  0  ,  1-a  >_  l-e‘ 
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In  a  similar  way  to  Corollary  3.7  we  obtain 


Corollary  3.9 

If  F  is  SIFR,  symmetric  about  the  origin,  1-a  >  1-e  *  and 

_  r-r'-l 

^  m  , 

q  >_  1-e  ,  then 


H  (r ' , r) 0  ,  (X) 

a ,  q  r  ,  r ,  n  — ' 


dF(x)  >  q  >  >  1-a 


-H  (r',r)e  .  (X) 

a,q  r  ,r,n  — 


It  can  be  seen  from  Corollary  3.9  that  we  generally  have  to  truncate  the 
sample  in  order  to  have  a  sharp  tolerance  limit.  For  example,  if  we  are 
computing  the  tolerance  limits  between  which  95%  of  the  symmetric  SIFR 
distribution  F  lies  with  99%  confidence  when  n  *  12  ,  we  will  need  to  set 
r'  *  2  and  r  -  11  in  order  that  the  limits  be  sharp. 

We  now  give  the  tolerance  limit  on  the  tail  probability,  1  -  F  ,  for 
values  to  the  right  of  the  median.  We  let 


H  (r'.r)  if  X^2(r-r'-l)  >_  -2(n-2)  log(l-q) 


C  (r\r) 
a,q 


~  if  X^2 (r-r'-l)  <  -2(n-2)  log(l-q)  . 


Theorem  3.10 

If  F  is  SIFR,  F (0)  -4,1-  F(x)  >_  F(-x)  for  x  >_  0  and 
0  <_  1-q  <_  4  ,  then 


PF  I1  -  F[Cl-«.2q-1(r,-r)5r'.r,„®]  i1'"} 
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Proof 


Let  G  be  the  Laplace  distribution,  then  from  Theorem  3.1 

f[h.  .(r',r)0  ,  (X)l  <*  gTh.  „  .(r'.rje,  (Y)l 

L  l-a,2q-l  r  ,r ,n  — J  -  (_  l-a,2q-lv  r',r,n-'J 

for  (n-2)H,  _  .(r',r)  <  1  and  <  q  <  1  .  It  follows  that  since 

1-a, /q-i  —  —  — 


we  have  proved  the  theorem  in  the  range 


(n-2)Hl-a,2q-l(r,’r)  -  1  * 


that  is 


Xl-a2(r'r,~1)  iL  -2(n~2)  log2 (1-q)  . 

Now  suppose  that  (n-2)H  „  . (r’,r)  >  1  •  By  noting  that 

1—Ot  |  fcVj*' J. 

F  r  — e  ,  (X)l  f^r^-e  ,  (Y)l 

L  n-2  r  ,r,n  —  J  -  L  n-2  r\r,n'-J 

iG[Va,2q-l<r''r)V,r,n®]  • 

the  proof  follows.  |  | 

2 

We  note  that  the  tolerance  limit  is  sharp  only  when  ^2(r-r'-l)  >_ 
-2(n-2)  log2(l-q)  .  The  following  corollary  gives  bounds  on  q  and  a  such 
that  this  condition  is  satisfied. 


Coro  1 Lary  3 . 11 


If  F  is  SIFR,  F(0)  =  ’i  ,  ] 
r-r'-l 


and  1-q  >  hi  e 


n-2 


,  then 


f(x)  *  F(-:  )  for  x  >  0  ,  1-n  ■*  1-e 
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\ 


PF{1-F[Hl-a,2q-l<r'-r)5r',r,n®]  i1^}  i  ^ 


I 


Proof 

By  Theorem  3.10  we  have  only  to  show  that 

XJ_a2(r-r'-l)  >_  -2(n-2)  log2(l-q)  . 

Let  K  denote  the  chi-squared  distribution  with  2(r-r'-l)  degrees  of 
freedom.  Since  K  is  IFR,  we  have  from  Barlow  and  Proschan  (1965)  that 

K[2(r-r'-l)]  <  1-e-1  . 


Now 


K[*i-«2(r-r,-i>]  •i-a 


therefore 


><L2(r-r'-i) 

2 ( r— r ’ — 1) 


>  1  . 


_  r'-r-l 

If  1-q  >_  %  e  n  ^  ,  we  have  proved  the  result.  |  | 

In  noting  that  F(0)  -  H  by  assumption,  we  see  from  Corollaries  3.7 
and  3.11  that  we  can  obtain  sharp  tolerance  limits  on  1  -  F  for  most  q  , 
and  a  values  of  interest. 

If  F  is  SIFR  and  symmetric  about  the  origin,  we  have  from  Theorem  3.10, 
Corollary  3.12 


If  F  is  SIFR  and  symmetric  about  the  origin,  then 
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C.  (r',r)6  ,  (X) 

l-a,q  r ,r,n  — 

p  \  J  dF(x)  1  q)  i  1-a  . 

-C.  (r ' ,r)0  ,  (X) 

l-a,q  r  ,r,n  — 


Constructing  a  proof  similar  to  Corollary  3.11,  we  obtain  Corollary  3.13 
for  the  symmetric  case. 


Corollary  3.13 

If  F  is  SIFR,  symmetric  about  the  origin,  1-a  >_  1-e  ^  and 

_  r-r’-l 

i  n-2  , 

q  1-e  then 


H.  (r',r)0  ,  (X) 

1-a , q  ’  r  , r , n  — 


dF(x)  q\  >_  1-a 


-H.  (r',r)0  ,  (X) 

l  1— ci , q  r  ,r,n  — 


We  see  from  Corollary  3.13  that  the  tolerance  limit  in  Corollary  3.12 
is  sharp  for  most  q  and  a  values  of  interest. 

Tolerance  limits  for  SDFR  distributions  although  probably  not  as  useful 
as  those  for  SIFR  distributions  are  nevertheless  quite  interesting.  Again, 
we  give  one  sided  tolerance  limits  when  F  is  skewed  and  two  sided  limits 
when  F  is  symmetric. 


Theorem  3.14 


If  F  is  SDFR,  F(0)  -  H  ,  1  -  F(x)  _>  F(-x)  for  x^0,0^q^»s 

and  (n-2)H  .  „  (r',r)  <  1  then 

a,l-2q  — 


M1-FK.l-2q,r'ir>5r\r.„<i)]  -  1-0 
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Proof 


Let  G  be  the  Laplace  distribution,  then  by  Corollary  3.2 

f[-H  .  ,  (r\r)0  ,  (X)l  g[-H  ,  „  (r\r)6  ,  (Y)l 
L  a.l-2q  r  , r , n  —  J  -  [_  ct,l-2qv  r\r,n-'J 


for  (n-2)H  ,  »  (r1 ,r)  <  1  .  By  noting  that 

a ,  l-2q  — 

Pc{GKl-2q(r’>r)5r\r.„©]  i’}  ■  l'“ 

we  have  proved  the  theorem.  | | 

In  a  similar  way  to  Corollary  3.7  we  see  that  if  1-a  ^  1-e  *  ,  then 
for  the  above  tolerance  limit  to  exist 


q  *5  exp 


X*2(r-r'-l)  | 

—(1^2)—  |  i  15  exp 


r-r'-l 

n-2 


>  h  e 


-1 


Thus  we  see  that  we  cannot  achieve  high  coverage  with  this  tolerance  limit. 

Now  in  the  case  that  F  is  symmetric  we  can  obtain  two  sided  tolerance  limits. 
From  Theorem  3.14  we  obtain  immediately 


Theorem  3.15 


If 


F  is  SDFR,  symmetric  about  the  origin,  and 


(n-2)H  (r'  ,r)  <  1  , 

a,q  — 


then 


H  (r\r)e.  (X) 

a,q  r  ,r,n  — 


dF(x)  >_  q\  l-o 


-H  (r ' ,r)6  ,  (X) 

a,q  r  ,r,n  — 


Again  in  a  similar  way  to  Corollary  3.7,  we  see  that  if  1-a  >_  1-e  ,  then 

for  the  above  tolerance  limit  to  exist  it  is  necessary  that 
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q  <_  1  -  exp 


I  (  r-V-A  ,  -1 

\  *  20.-2) - )  i1  -  i1'6 


and  so  high  coverage  Is  not  possible. 


Theorem  3.16 


If  F  is  SDFR,  F(0)  -4,1-  F(x)  >_  F(-x)  for  x  >  0  ,  >s  <  q  <  1  , 


and  mH 


fll-a,2q-l(r'*r)  -  1  »  then 

Pf{X  -  F[Hl-a,2q-l<c'-r>V,r,n®]  i  ^  }  1 


Proof 


Let  G  be  the  Laplace  distribution,  then  by  Corollary  3.2, 

F[Hl-<,.2q-l‘''®V,r.n®]  I*  C[-Hl-a.2<,-l(t'-r>V,r,n®] 

for  mH.^  2q-l^r',r^  —  ^  noting  that 

Fc{G[Kl-».2q-l(r’>r)5r',r,n®]  i’}  ' 


we  have  proved  the  theorem.  |  | 

If  1-a  >_  1-e  *  ,  then  by  similar  reasoning  to  that  used  in  Corollary  3.11, 
we  see  that  for  the  above  tolerance  limit  to  exist 


/  *L2<r-r'-l)\  .  /  r-r'-l\ 

l-q  2.  ^  exP  ^ - 2S - )  i  >s  exp  I - ^ - 1  . 

It  will  thus  be  necessary  to  truncate  the  sample  rather  severely  in  order 
that  the  tolerance  limit  on  the  tail  probability  exists. 

If  F  is  SDFR  and  symmetric  about  the  origin,  then  from  Theorem  3.16 


we  have , 
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In  the  same  way  as  for  Theorem  3.16,  the  sample  will  have  to  be  truncated 
rather  severely  In  order  that  the  tolerance  limit  exists. 
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3.2  Tolerance  Limits  Based  on  the  Standard  Deviation 

Frequently  sampling  plans  are  based  on  the  assumption  that  the  underlying 
distribution  Is  normal,  which  gives  rise  to  tolerance  limits  of  the  form 
x  +  A 1 s  where  x  and  s  are  the  sample  mean  and  standard  deviation 
respectively.  If  the  mean  p  is  known,  which  we  will  assume,  the  tolerance 
limits  are  of  the  form  p  +  *£8  .  Since  we  are  generally  not  sure  that  the 
assumption  of  a  normal  distribution  is  correct,  it  is  interesting  to  see  how 
robust  the  normal  tolerance  limits  are  against  s-ordered  alternatives.  The 
following  theorem  throws  some  light  on  this. 

Theorem  3.18 

If  F  <  G  and  F  ,  G  are  symmetric  about  the  origin,  then 
8 


for  all  a  e  [0,1]  and  A  <_  1  ,  (A  >  /n)  ,  where 


Proof 

Define  G2  by 

G2(x)  =  P(Y2  i  x) 

2  .  We  see  that  “  [G  *F(y)]2  and  by 


and  similarly  for  F 
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differentiation  we  find  that  ^^(x)  is  convex  if  G  ^F(x)  is  convex 
for  x  >  0  .  From  Lemma  3.1  of  Barlow  and  Proschan  (1966a) 


<(>) 


for  X2  <_  1  (X2  >_  n)  .  Since  G^F^X2)  Y2 


_  ,,2  2,  st  ,st.  _  ,,2  2. 
^2 ( X  s^)  <_  (>_  )  G2 ( X  Sy) 


The  theorem  now  follows  from  the  observation  that  F2(x)  *  F(i'/x)  -  F(-/x)  .  |  | 
Applying  the  strong  law  of  large  numbers  we  obtain 


Gorollary  3.19 

If  F  <  G  and  F  ,  G  are  symmetric  about  the  origin,  then 
s 


Xo. 


Xo, 


dG(x) 


for  X  <  1 


-Xo 


-Xo, 


and  o  is  the  variance  of  X  • 
X 


This 


is  the  symmetric  analogue  of  F(0^)  <_  G(02)  ,  when  F  <  G  and 


0^  ,  ©2  are  the  means  of  F  and  G  respectively.  [See  Barlow  and 
Marshall  (1964),  Theorem  7.1.] 

As  an  example  of  Theorem  3.18  we  let  G  be  the  normal  distribution  with 
known  mean  but  unknown  variance.  If  we  consider  the  tolerance  limits  p  +  Xs 
for  G  which  contain  99%  of  the  population  with  99%  probability  we  see  that 
if  the  number  of  observations  from  G  is  n  then  X  >_  /rT  for  n  18  . 
Hence  we  can  assert  that  if  F  <  G  ,  F  symmetric  and  G  the  normal  and  both 
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with  known  means,  then  for  a  sample  size  less  than  18,  the  99%  tolerance  limit 
with  99%  confidence  based  on  the  normal  is  conservative  for  F  . 
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3.3  Confidence  Bound  on  the  Variance  of  a  Symmetric  SIFR  Distribution 

In  order  to  calculate  a  confidence  bound  on  the  variance  of  a  symmetric 
SIFR  distribution  we  need  to  develop  bounds  on  the  distribution  in  terms  of 
the  mean,  which  we  will  assume  to  be  zero,  and  the  variance.  We  note  that 
two  symmetric  SIFR  distributions  with  the  same  mean  and  variance  must  cross 
at  least  three  times.  Bounds  will  only  be  given  to  the  right  of  the  mean, 
bounds  to  the  left  following  by  the  symmetry  of  F  . 


Lemma  3.20 


If  F  is  SIFR  and  symmetric  about  the  origin  with  variance  o  ,  then 


F(x)  > 


0  <  x  <  o 


i  i  -bx 

,1  -  Hi  e  x  >  a 


where  b  is  a  function  of  x  and  o  and  is  given  by  the  solution  to 


2 

o 


and  the  bound  is  sharp. 


Proof 


The  bound  for  0  <_  x  <_  o  is  obvious  and  is  attained  by  the  symmetric 
distribution  with  ro  iss  ^  at  -o  and  +o  .  For  x  >_  o  we  consider  the 
distribution  G^,(x)  which  is  symmetric  about  the  origin  and  for  positive 
arguments  is  given  by 


1  -  h  e~bx  0  <  x  <  T 


Gt(x) 


x  >  T 
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where  b  is  determined  by 


2 


oo 

/ 


x2dGT(x) 


2 

o 


We  note  that  since  F  must  cross  G  at  least  once  in  (0,°°)  ,  and  G_1F(x) 
is  convex  for  x  >_  0  that  F(x)  >_  GT->x(x)  ,  x  ^  o  .  j  | 

Lemma  3.21 

If  F  is  SIFR  and  symmetric  about  the  origin  with  variance  o2  ,  then 


sup  e 

0  <  A  <  o 


-a(x-A) 


0  <  x  <  o 


F(x)  < 


X  >  0 


where  a  is  given  by  the  solution  to 


2/a2-o2\ 

a  i  -  I  +  aA  +  1  =  0 


Proof 


Consider  the  distribution  G^Cx)  ,  symmetric  about  the  origin,  and  which 
for  positive  arguments  is  given  by 


h 


ga(x)  = 


1  -  L  e-a(x-A) 


0  <  x  <  A 


x  >  A 


and  where  a  is  given  by 
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2 


/ 


x  dG . (x) 

A 


2 

o 


and  A  lies  in  the  range  0  <_  A  <_  o  .  Since  (x)  is  SIFR  with  mean  zero 
2 

and  variance  o  ,  it  must  cross  F  at  least  once  in  the  interval  (0,°°)  . 
Let  G^  cross  F  from  above  at  x  *  u(A)  .  If  u(A)  does  not  exist  set 
u(A)  *  °°  .  Let  us  consider  the  following  two  cases: 

Case  (i)  u(0)  >_  a 

Since  G  ^F(x)  is  convex  for  x  >_  0  ,  we  see  that 

F(x)  <_  Gq(x)  0  1  x  1  0  • 


Case  (ii)  u(0)  <  o 

Since  G  ^F(x)  is  convex  for  x  ^  0  ,  we  see  that  u(A)  t  A  ,  u(o)  *  °° 
and  u(A)  is  continuous  in  A  on  the  positive  support  of  F  .  Therefore 

_ .  .  -a(u-A)  _ 

F(x)  sup  e  0  <_  x  o 

0<_A<_o 

The  bound  for  x  >  o  is  obvious  and  is  achieved  by  the  distribution  with 
mass  4  at  x  ■  -o  and  x  =  +o  .  |  [ 

We  are  now  in  a  position  to  derive  a  lower  confidence  bound  on  the 
variance  of  a  symmetric  SIFR  distribution. 

Theorem  3.22 


If  F  is  SIFR  and  symmetric  about  the  origin,  then 
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From  Lemma  3.20  we  have  the  bound 


t  <_  o 

t  >  a  , 

and  b  is  given  by  the  solution  to  (3.1).  If  G  is  the  Laplace  distribution, 
then  since  F  is  SIFR 


Since 


G  Mr  0r,  r  n(Y)\  <  1 

I  yn-2  r  ,r,n  J 


xl-a2 (r-r ' -1) 


4  e  2(n“2) 


we  have 


since  o  I 


b  .  Hence  we  have  the  result.  | 
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IV.  PROPERTIES  OF  THE  ORDER  STATISTICS 


Marshall,  Olkin  and  Proschan  (1965)  and  Barlow  and  Proschan  (1966a)  have 
developed  Inequalities  for  order  statistics  arising  from  distributions  F 
and  G  in  the  case  that  G  ^"F  is  starshaped  on  the  support  of  F  and 
G(G  )  »  F(0  )  *  0  .  Van  Zwet  (196^)  has  extensively  treated  inequalities  for 
the  expectations  of  order  statistics  arising  from  e-ordered  and  symmetric 
s-ordered  distributions.  These  inequalities  are  not  only  interesting  but  are 
useful  in  developing  bourJ  .  and  giving  insight  into  the  nature  of  convex  and 
concave-convex  transformations.  We  shall  develop  inequalities  for  the  order 
statistics  and  for  power  combinations  of  the  random  variables  in  the  case  of 
two  symmetric  r-ordered  distributions,  except  for  one  inequality  where  we 
require  the  stronger  s-ordering.  The  inequalities  reflect  the  fact  that  an 
antisymmetric  starshaped  transformation  of  a  random  variable  shifts  mass  to 
the  tails. 


Barlow  and  Proschan  (1966)  showed  that  if  F(0)  *  G(0)  -  0  and  — — 

is  nondecreasing  in  x  ^  0  ,  then  the  ratio  of  order  statistics  EY^  ^/EX^  n 

is  also  increasing  in  i  »  1,2,  ...,  n  for  all  n  ;  i.e.,  r-ordering  on  the 

positive  axis  is  preserved  by  the  expected  values  of  the  order  statistics. 

Van  Zwet  (1966)  showed  that  if  F  <  G  then  JEY , , ,  -  EY .  ]/[EXJt1  -  EX.  ] 

c  i+l,n  i,nJ  i+l,n  i,nJ 

is  nondecreasing  in  i  '  1,2,  ...,  n  for  all  n  .  He  also  proved  that  if 
F  <  G  ,  F  and  G  symmetric  about  the  origin,  then  the  expected  values  of  the 
order  statistics  are  similarly  s-ordered.  Independently  of  van  Zwet  we 
proved  a  related  result,  namely  that  for  r-ordered  symmetric  distributions  the 
expected  values  of  the  order  ?tistics  preset ve  the  ordering. 

We  shall  need  the  concept  of  total  positivity.  A  function  K(x,y)  of 
two  real  variables  x  c  X  ,  y  e  Y  ,  where  X  and  Y  are  ordered  sets,  is 
said  to  be  totally  positive  of  order  r  (TP^)  if  for  all  1  <  n  <  r  , 
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x.  <  x„  <  . . .  <  x  and  y,  <  y„  <  ...  <  y  ,  where  each  x.  t  X  ,  e  Y 
1  —  l  —  —  m  1—  2—  —  m  l  1 

we  have  the  determinantal  inequalities 

m 

K(x  ,y  )  >0. 

2  i» j*l 

For  further  treatment  of  TP  functions  and  their  properties,  see  Karlin  (1964). 
It  is  convenient  to  define 

K(1,n’X)  '  ('i~l)f (n-i)T  Fi_1(x)  • 

We  will  need  the  following  well-known  result: 

Lemma  4 . 1 

i\  (i,n,x)  is  TP^  in  i  *  1,2,  ...  and  -<*>  <  x  <  00  . 

?r<^of 

Since 


is  a  generalized  Vandermonde  matrix,  we  know  (see  Gantmacher  (1939),  p.  118) 
that  it  is  totally  positive.  The  lemma  follows.  |  | 


t  n  >  i 
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(ii) 

Proof 

Let 


EX  ...  EX, 
n-i+1 tn  i,n 

EY  ...  “  EY, 

n-i+1, n  i,n 


h(i) 


00 

/ 


(x  -  cG_1F(x))K(i,n,x)dF(x) 


EY 


‘•fo ' c) 


Since  F  is  symmetric,  we  have  that 


h(i) 


-h(n-i+l) 


for 


i.e.,  h(i)  is  antisymmetric  about 


n+11 


Since  K(i,n,y)  is  TP  for 


i  -  1,2,  x  <  00  ,  and  (x  -  cG  ^F(x))  changes  sign  at  most  three 

times  for  c  ^  0  ,  we  have  by  the  variation  diminishing  pruperty  of  TP^ 
functions  that  h(i)  must  change  sign  at  most  three  times.  If  h(i)  does 
change  sign  three  times  then  the  order  of  the  signs  must  be  the  same  as  for 


(x  -  cG  *F(x))  ;  viz,  +  -  +  - 

EX  r  - 

we  see  that  —  >n  l  i  >  NJ. 

EY,  —  12 

i,n  L  ■ 


Since  h(i)  is  antisymmetric  about 
+  1  ,  proving  (i). 


n+1 


(ii)  may  be  most  readily  proved  by  a  geometrical  argument. 
By  conditioning  on  the  (n+l)st  observation,  we  see  that 


ex,  -  (i  -  -rrlEx,  +  -4r  ex_.  ..  , 

i,n  \  n+ly  i,n+l  n+1  i+l,n+l 


and  hence  for  n  <  2i  we  obtain  the  following  diagram  from  which  we  observe 
EY 

—  J~n  i  n  for  n  .<  2i  . 

EX, 

i,n 
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Similarly  we  can  show  that 


EY 


EX 


I  n  for  n  ^  2i 

i,n 


We  are  now  in  a  position  to  obtain  bounds  on  EX,  .  Note  that 

i,n 


EX 


EX 


l.i 


EX 


i.n 


7  +  I." 


EY.  -  EY.  -  EY 
i,i  i,n 


f  +  l,n 


for  i 


2.  +  1  ,  and  where  is  the  smallest  integer  larger  than  or 

equal  to  -j  .  The  first  inequality  follows  from  Theorem  4.2  (ii),  and  the 
second  from  Theorem  4.2  (i).  If 


e 


co 


x| dF(x) 


t 


then 


0EY 

'ey 


ixH 


<  EX 

i.n 


6EY 


LlD _ 

+  l,n 


for  i  >_  ~  +  1 


Note  that  no  nontrivial  upper  bound  exists  when  n  is  odd. 
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Similary  for 


-6EY 


EYr- 


i.n 


-0EY 


n 

2  »  n 


<  EX 


“  i,n  -  EY 


i,n 


l,n-i+l 


Using  Theorem  4.2,  inequalities  for  linear  combinations  of  order 
strtistlcs  from  r-ordered  distributions  may  be  derived  in  much  the  same  way 
as  Barlow  and  Proachan  (1966a)  do  for  star-shaped  ordering  (i.e., 

F(0)  -  G(0)  -  0  and  G  1F  starshaped  on  the  support  of  F) .  As  these 
inequalities  parallel  those  of  Barlow  and  Proschan  (1966a)  we  shall  omit  them. 

Van  Zwet  has  obtained  necessary  and  sufficient  conditions  on 
al*a2»  •••»  an  such  that 

for  F  <  G  .*  We  derive  sufficient  conditions  on  a.  . . a  such  that 

-1-4  n 

the  above  inequality  is  true  when  F  <  G  .  We  strongly  believe  that  these 
conditions  are  also  necessary  but  have  not  been  able  to  show  this. 

Theorem  4.3 

If  F  <  G  ,  G  symmetric,  F(0)  -  G(0)  -  h  and  1  -  F(x)  F(-x)  ,  x  >  0 

then 


(4.1) 


'(f  °iEXl.n)  i  °([  *!«!,„) 


when 


■ft 

Personal  communication. 
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(4.2)  0  <  Y  Aj+1[yJa-y)n'j  +  yn"j(i-y)J]-^i  -  1 


and 


(4.3) 


n-l  /  \ 

J-0  (j) 


Vi  yn'J(1'y)3  -*>  'i  0 


for  all  h  <  y  <  1 


Proof 

We  prove  (4.1)  for  the  unbracketted  inequalities  first. 

Let 


h(x)  “  l  ai  (l-l)?in~i)T  G1"1(x)(l-G(x))n"iG’(x)  . 


.1-1, 


Now  by  letting  4,  =  -F  ^G  the  theorem  can  be  re-expressed  as  follows.  For 
<f>  concave-convex  about  the  origin  and  <>(x)  <_  -$(-x)  for  x  _>  0  and  G 
is  symmetric  about  the  origin  then 


(4.4) 


/  ?  \  ? 

4>f  /  xh(x)dxl_^  /  <}>(x)h(x)dx 

\— oo  /  —00 


if  (4.2)  and  (4.3)  are  satisfied  for  all  *5  <  y  <  1  .  By  an  argument  similar 
to  that  used  in  Theorem  2.1  we  can  see  that  (4.1)  is  true  if  and  only  if  (4.4) 
holds  for  all  x^  xQ  _>  0  ,  when  <j>  is  a  double  angle  function:  i.e.. 


x  +  xQ  for  x  -Xy 

4>(x)  0  for  -xy  <  x  <  x^ 

x  -  x^  for  x  _>  x^ 
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For  $  a  double  angle  function,  (4.4)  becomes 


T 


(x+Xg)h(x)dx  + 


00 

/ 


)h(x)dx 


(4.5) 


Imax 

0,  J  xh(x)dx  -  x  J 

(mini 

(  -L  / 

f  \ 

0,  1  xh(x)dx  +  x0J 

\  ' 

l  / 

+00 

for  J*  xh(x)dx  >_  0 


oo 

/ 


Let  p (y)  -  l  ai  Yn-i)"(i-l) '  y^1^)"  1  and  substituting  y  -  G  (x) 
in  (4.5)  we  obtain 


i,,  .n-i 


i-yo  i 

f  (G_1(y)  -  G"1(l-y0))p(y)dy  +  j  (G_1(y)  -  G"1(y1))p(y)dy 

n  J. 


(4.6)  I  max 


f!‘“ 

fiv 


(y)p(y)dy  -  g' 


1(y1)|  for  J ' G  1(y)p(y)dy  _>  0 


,-l. 


(y)p(y)dy  -  G  (1— y0)|  for  /  G  J'(y)p(y)dy  <  0  , 


fo\ 


and  (4.1)  is  true  if  and  only  if  (4.6)  is  true  for  all  1  ?  y^  yg  4 

and  all  G  ^  strictly  increasing  on  (0,1)  and  antisymmetric  about  h  . 
We  shall  now  show  that  (4.2)  and  (4.3)  together  imply  (4.6). 

From  the  well-known  equality 


'(n-i)i(i-i)!  yV'W'V  ■  l  (")  .J(l-n)n-3 

n  “ 
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Since  any  strictly  increasing  function  defined  on  (0,1)  and  antisymmetric 
about  may  be  approximated  arbitrarily  closely  from  below  (above)  for 
x  ^  h  (x  <  !i)  by  a  positive  multiple  of  double  step  functions,  and  since 
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for  all  h  <_  y  1  . 

Van  Zwet  (1964)  proves  that  if  F  <  G  and  F  ,  G  symmetric  about  the 
origin,  then 


[e  !  X ;  “]  [E|Yla] 


for  those  values  of  a 

such  that 

Ely;3 

exists  and 

...  2k 

_,,2k 

(4.10) 

EX 

■  ■  ■  < 

9  k  - 

EY 

for  k  *  1,2, 

(EX  ) 

(EY") 

2k 

for  those  values  of  k  for  which  EY“  exists.  We  will  prove  a  stronger 
result;  namely  that  given  F  <  G  then  the  inequalities  (4.9)  and  (4.10)  hold 
stochastically  for  the  usual  estimates  of  the  expectations,  and  hence  by  the 
strong  law  of  large  numbers,  for  the  expectations  themselves. 

We  need  to  introduce  the  concept  of  majorization ,  and  one  of  the  theorems 
applying  this  concept.  For  a  fuller  treatment  see  Hardy  Littlewood  and  Polya 


(1959)  and  Ostrowski  (1952). 
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Definition 

(i)  A  sequence  a  ■  (a^,  a^)  is  said  to  majorize  a  sequence 

b  ■  (b, ,  ....  b  )  (written  a  >*  b)  if  a,  >  . . .  >  a  ,  b,  >  . . .  >  b  ,  and 

r  r  n  n 

l  a  _>  l  b  for  r  -  1,  ...,  n-1  ,  while  £  a  =  £  b  . 

1  1  1  11 

Theorem  4 . 7  (Hardy,  Littlewood  and  Pdlya) 

If  <$>  is  convex  on  the  interval  I  and  where  x^,  ; 

yl*  *•••  yn  beion8  to  1  »  then 


l  1 1 

i  i 


Theorem  4.8 


If  F  <  G  ,  F  and  G  symmetric  about  the  origin,  then 


1WT  J"4 

*n  "lb  —  rn  "l  1 

I  |x±|a  l  |Y±|a 

Li  J  Li  J 


and  if  E  Y  exists  then 


>11  f _  <  X  .  I 

b  —  b 

( E  |  X |  a)  (E  |Yl  a) 


Proof 


0  <  a  <  b 


Raise  to  the  a  power  the  absolute  value  of  the  observations  from  F  , 


and  order  so  that 


|X|*_  |X|*  >  ...  >  |X|- 
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i 


I 


l 


% 


Now  if  F  (x)  *  P  (  |  X 1 3  <  x)  ,  a  >  0  we  see  that 

3  r 

since  Q  --  t  x  >  0  we  have 

x  — 


G'1?  (x3) 
a  a 


[G_1F(x)]a  and 


F  (x) 
a _ 


x 


t  x  0 


If  follows  from  a  theorem  in  Marshall,  Olkin  and  Proschan  (1966)  [cf.  Barlow 
and  Proschan  (1966),  Theorem  3.12]  that 


for  k  =  1,2,  ....  n  . 

Q 

Now  from  Theorem  4.7,  by  considering  the  convex  function  <f> (x )  *  x  , 
x  ^  0  ,  c  1  we  obtain  the  stochastic  inequality 


I  !*i 


|3C 


n 

V 

St  1 

< 


,ac 


(!  "  (i  lv‘'f 


Letting  ac  *  b  we  obtain 


(I  "■i>)~  JM 

a  ;  a ».) 


Now  if  E | Y | b  exists,  then  E | Y [ a  exists  and  by  a  limiting  argument 
we  can  see  that  E I X j b  exists.  (ii)  is  then  true  by  the  strong  law  of  large 
numbers. 


5 


Corollary  4.9 

If  F  <  G  ,  F  and  G  symmetric  about  the  origin,  then 


(1) 


l  ^  t  I  Y^rk 

m'¥! 


r  -  1,2,  ... 
k  -  1,2,  ... 


2rk 

(11)  If  EY  exists  then 


EX' 


2rk 


EY 


2rk 


2r  k  2r  k 

(EX  r)  (EY  ) 


Corollary  4,10 


2rk 


If  F  <  G  ,  F  and  G  symmetric  about  the  origin  and  EY^  exists, 


EX 


2rk 


EY 


2rk 


2r  k  2r  k 

(EX^r)  (EYp 


for  k  ■  1,2,  . . . 


Proof 

The  proof  follows  in  the  same  way  as  for  Corollary  4.9  and  by  the 
observation  that  G(£)F(i)^x^  “  G  *F(x)  »  where  F^  (x)  -  P(X^  <_  x)  . 


then 
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